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Abstract 

We have measured the thermal properties of suspended membranes from 10 K to 300 K for two 
amplitudes of internal stress (about 0.1 GPa and 1 GPa) and for two different thicknesses (50 nm 
and 100 nm). The use of the original Sca-Volklein method has allowed the extraction of both the 
specific heat and the thermal conductivity of each SiN membrane over a wide temperature range. 
The mechanical properties of the same substrates have been measured at helium temperatures using 
nanomechanical techniques. Our measurements show that the thermal transport in freestanding 
SiN membranes is not affected by the presence of internal stress. Gonsistently, mechanical dissipa¬ 
tion is also unaffected even though Qs increase with increasing tensile stress. We thus demonstrate 
that the theory developed by Wu and Yu [Phys. Rev. B, 84 , 174109 (2011)] does not apply to 
this amorphous material in this stress range. On the other hand, our results can be viewed as a 
natural consequence of the ’’dissipation dilution” argument [Y. L. Huang and P. R. Saulson, Rev. 
Sci. Instrum. 69 , 544 (1998)] which has been introduced in the context of mechanical damping. 
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I. INTRODUCTION 


Silicon nitride (SiN) thin films are widely used to thermally isolate sensitive thermal 
detectors, for etch masking as well as layers for micro-electromechanical systems^. Indeed, 
outstanding mechanical properties including very high quality factors can be reached 
in optimized SiN material. Depending on deposition parameters, SiN films can experience 
very large residual (biaxial) stress during deposition. It is thus of prime importance to 
understand the role of the internal stress not only on the mechanical properties but also 
on the other physical characteristics of SiN films including optical, thermal and electrical 
properties. Silicon nitride has a specific place due to its amorphous nature and the study 
of stress in that compound is also an issue for the fundamental understanding of its role in 
the physics of glasses^. 

Using the stress to tune the thermal properties of nanomaterials is one of the possible 
ways to design future thermal components (thermal rectifier, thermal diode, thermal switch^ 
...). This has been proposed for monocrystalline silicon®’®, as strain in silicon is currently 
used to enhance electron mobility in transistors^. Since the debate on the origin of mechan¬ 
ical dissipation in strained glasses like SiN^’® the question of the effect of stress, whatever 
its origin (internal or external) on the thermal properties has been raised and theoretically 
addressed for the case of silicon nitride®. Indeed it is well known that stoichiometric silicon 
nitride (Si 3 N 4 ) prepared by low pressure chemical vapor deposition (LPCVD) contains a sig¬ 
nificant internal tensile stress (up to about 1 GPa) as compared to regular non-stoichiometric 
SiN that has a very low internal stress (below 0.2 GPa). In an attempt to explain the very 
high mechanical Qs, Wu and Yu® proposed a model where the internal losses in the material 
are sensitive to the stress state. Their calculations based on this hypothesis predict that 
the thermal conductivity of SiN may be strongly enhanced by the presence of stress. On 
the other hand, systematic mechanical measurements on high stress SiN substrates explain 
the Qs through the so called ’’dissipation dilution” modefi®: mechanical energy is stored 
through the tensioning of the substrate while the dissipation is unaffected®’^^’^®. However, 
no experiments to date compared directly similar devices made of different SiN materi¬ 
als. Furthermore, thermal properties of low stress SiN have been widely measured over a 
broad temperature range^®^^^ for different kinds of thin films and nanomaterials, but very 
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few experimental studies deal with the influence of stress on the thermal transport at the 
nanoscale 

In order to study the potential effect of internal stress on the thermal properties of silicon 
nitride, the thermal conductivity and the specihc heat have been measured as a function of 
temperature for high stress (HS) and low stress (LS) SiN membranes having a thickness of 
50 nm and 100 nm. These measurements are performed using the Sca-Volklein method^®^^^ 
as described in previous papers. This technique allows the measurement of both thermal 
conductivity and specihc heat of a given membrane within the same experiment over a 
broad temperature range. The mechanical dissipation and stress amplitude of both HS and 
LS substrates are also measured at cryogenic temperatures by means of nanomechanical 
resonators^^. We show experimentally that thermal conduction is essentially independent 
of the stress stored in this material. This is inconsistent with the hypothesis underlying 
the model of Ref.® for SiN, and corroborates the ’’dissipation dilution” explanation for high 
mechanical Qs. 

II. SAMPLES AND EXPERIMENTAL METHODS 

The thermal properties of two types of SiN membranes have been measured: high stress 
stoichiometric Si 3 N 4 and low stress SiN deposited by LPCVD. The amorphous stoichiometric 
high stress (HS) Si 3 N 4 as well as low stress (LS) SiN were grown on both sides of a silicon 
substrate. The membranes were then patterned on the rear side by laser photolithography. 
After removing the silicon nitride by SFe Reactive Ion Etching, the silicon substrate on the 
rear side was etched in KOH, as described in Fig. 1. The hnal result is a rectangular SiN 
membrane obtained on the front side. 

Before the thermal study, a mechanical measurement was performed to quantify the 
stress present after releasing the membranes. A suspended silicon nitride beam with 100 nm 
thickness, 250 nm width and 15 pm length fabricated using e-beam lithography from the 
same substrate was placed in a magnetic held (see Fig. 1 for fabrication details). A sinusoidal 
driving current within a 30 nm thin deposited A1 layer is used to generate the Lorentz force 
causing the beam’s out-of-plane oscillation^^. This measurement is performed in a vacuum 
of about 10“® mbar at helium temperatures. The magnetic hux cut by the beam oscillation 
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FIG. 1. (colour online) Micro and nanofabrication processes of both suspended structures studied 
in this work, a) fabrication process for the membrane 1), 2) The patterns of the membranes are 
created by photolithography. The non protected SiN is removed by SFg RIE etching. 3) The silicon 
is anisotropically etched in a KOH solution. 4) The thermometers are obtained by a lift-off process; 
the area is patterned by photolithography. 5) NbN (70 nm) is deposited by reactive sputtering. 
6) The resist and NbN layer is removed using a wet procedure, b) fabrication process for the 
nanowires 1) The patterns of the nanowires are created by e-beam lithography. 2) evaporation of 
A1 layer (30 nm) 3), 4) The non-protected SiN is removed by SFg RIE. 5) The silicon is isotropically 
etched by gazeous XeF 2 ecthing. 


generates a voltage which is measured using a lock-in ampliher^’^^. Typical resonance curves 
for the first flexure and their respective fits are shown in Fig. 2. 

The expression for the mode resonance frequency of a stressed doubly-clamped beam 
is given by^®: 




( 1 ) 


Eq. 1 is used to calculate the stress a within the beam, with n the mode number, p the 
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FIG. 2. (colour online) Measurement of the first flexure resonance of two suspended SIN beams 
with different in-built stresses, the red points are the in-phase signal and the black points the 
out-of-phase signal (both 15 /rm long, 250 nm wide, 100 nm thick, see SEM picture in inset), 
(a) low-stress and (b) high-stress resonance lines obtained in the linear regime. The lines are 
Lorentzian fits, with full width at half height of 650 ± 50 Hz (HS, Q ~ 25 000) and 500 ± 50 Hz 
(LS, Q ~ 14000). We extract from the resonance frequencies the stress values of 0.85 ± 0.08 GPa 
(HS) and 0.12 ± 0.05 GPa (LS). Data taken at 4.2 K in vacuum. 


silicon nitride density (3 g/cm^) and h the beam length^^. We find a stress value of about 
0.85 GPa for the HS silicon nitride, which conhrms that the membrane is still stressed after 
releasing and 0.12 GPa for the LS silicon nitride. Note that both values fall in the high- 
stress limit of beam theory, validating the use of Eq. (1). At the same time, the linewidth 
of the resonances in Fig. 2 (measuring mechanical dissipation) are almost equal for HS and 
LS beams. Even though the stress is modihed by almost an order of magnitude, mechanical 
friction around 4 Kelvin seems insensitive. Many devices varying shapes and stress have 
been measured, studying flexural modes up to n = 9 and conhrming these Endings^®. The 
temperature dependence of mechanical properties shall be discussed elsewhere^®. 

These mechanical measurements seem to be in contradiction with the claims of Ref.^ 
(in which very different types of devices were compared), but the raw data do agree: we 
End indeed a higher Q for a higher stress^’®. On the other hand, our results are a natural 
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consequence of the ’’dissipation dilution” modeh®^^^ in which high Qs are reached because 
of the flexural energy stored in the tensioning of the structures: while in Refsd^’^^ only HS 
devices are measured, we show here with a comparison between geometrically identical HS 
and LS structures that there is no internal effect of stress on mechanical dissipation. 

Thermal experiments are conducted on the very same materials. As mentioned above, we 
have chosen in this study the most appropriate method to measure the thermal conductivity 
of large aspect ratio suspended membranes: the Sca-Volklein method^®^^^. The principle of 
the method consists in creating a sinusoidal Joule heating generated by an AC electrical 
current at frequency u across a transducer centered along the long axis of a rectangular 
membrane. The center of the membrane is thermally isolated from the frame and hence its 
temperature is free to increase. 

The temperature oscillation (^100 mK) of the membrane is at 2u and is directly related 
to its thermal properties by the amplitude and the frequency dependence of the aforemen¬ 
tioned temperature oscillation. Since the resistance of the thermometer can be considered 
as linearly dependent on temperature over that small temperature oscillation, the voltage 
V = R[T{2u)] X I{u}) will have an ohmic component at u and a thermal component at 3u. 
By measuring the voltage appearing across the transducer as function of the frequency, 
it is possible to deduce the thermal conductivity and the specihc heat of the membrane^^. 
The membranes measured in this study are 300 pm wide and 1.5 mm long. 

The transducer of 5 pm width and 1.5 mm length is made out of NbN whose resistance 
is strongly temperature dependent. It serves as a thermometer and a heater at the same 
time30-32 Pqj, present work, the thermometer has been designed for the 10 K to 320 K 
temperature range. Typically, the resistance of the thermometer is about 100 kOhm at 
room temperature with a temperature coefficient of resistance (TCR) a = of 10“^ 
at 300 K and of 0.1 at 4 K. 

Since the lea voltage is 3 to 4 orders of magnitude higher than the 3a; voltage, a specihc 
Wheatstone bridge is used to reduce the la; component and perform thermal measurements 
(see Fig. 3). The bridge consists of the measured sample with a resistance Re, which is the 
NbN thermometer on the SiN membrane, the reference thermometer Rref deposited on the 
bulk region of the chip which has the same geometry and deposited in the same run as the 
transducer on the membrane, an adjustable resistor R^,, and an equivalent nonadjustable 
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FIG. 3. (colour online) Experimental set-up based on the Wheatstone bridge configuration, the 
yellow membrane sample is on the bottom right and the reference thermometer is on the left; the 
blue area of the sketch corresponds to the temperature regulated part of the Wheatstone bridge. 


resistance Ri =50 kOhm. 

The general expression of the measnred 3n; ontpnt Wheatstone bridge voltage can be 
given 




V^^^aReRi \/\R 


2u) 


2 (Re + Rl 


( 2 ) 


with a the TCR, Vac the lu input Wheatstone bridge voltage and |AT 2 a;| the amplitude 
of the temperature oscillation at 2oj of the membrane due to the sinusoidal nature of heating. 
By solving the partial differential equation of the heat flux across the membrane, eq. 3 


gives the relation between the thermal properties, the dimensions of the membrane and 

T/ 19,20. 






rms\3 JD jd 2 
ac ) 


4Tr2\ 11/2 


AKp (Re + Ri)' [1 + a;2 (dr^ + |^ + ^)] 


( 3 ) 


with Kp = ^ the thermal conductance and C = cSi the heat capacity of the measured 
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membrane, t = S- the thermalization time of the membrane to the heat bath, D = — 
the thermal diffusivity, i half the width of the membrane and S the section of the mem¬ 
brane (perpendicular to the heat flow). By measuring the voltage as a function of 
frequency both k (in-plane thermal conductivity) and c (specihc heat) of the membrane can 
be extracted^^. 



FIG. 4. (colour online) Thermal conductivity measurement of 50 nm and 100 nm thick membranes 
for both SiN low stress and high stress. The 100 nm curves of low stress and high stress show 
nearly no difference. 


III. EXPERIMENTAL RESULTS 

The thermal conductivity of the four different membranes (50 and 100 nm, low stress and 
high stress) has been measured versus temperature from 10 K to 300 K. The experimental 
data of thermal conductivity are presented in the Fig. 4. As expected for amorphous mate¬ 
rials the thermal conductivity of all membranes is continuously increasing with temperature 
as observed by Queen and Heilman^®. The general trend of the temperature variation of 
thermal conductivity of all different SiN membranes (LS and HS) are very similar. Only the 
50 nm LS membrane seems to have a slightly lower thermal conductivity at room tempera¬ 
ture with a value approaching 2.5 W.m“^.K“^, instead of 3 W.m“^.K“^ for the others. 

In all cases, values of the thermal conductivity at room temperature are approximately 



Reference 

Deposition 

Stoichiometry 

Stress 

k at 300 K 

c at 300 K 

Sample 

37 

LPCVD 

Sio.66No.34 

not measured 

3.2 

0.7 

free st. in- plane 

36 

LPCVD 

SiiNi.i 

not measured 

8-10 

not measured 

free st. out of plane 

33 

LPCVD 

Si rich 

low stress 

4.5 

0.5 

free st. in-plane 

34 

LPCVD 

Si rich 

low stress 

3.5 

not measured 

free st. in-plane 

35 

PECVD-APCVD 

SiiNi.i 

not measured 

0.3 

not measured 

out of plane 

38 

LPCVD 

not measured 

high stress 

3.2 

not measured 

free st. in-plane 

18 

LPCVD 

SiiNi.i 

from 0 to 2.4% 

2.7 (LS) to 0.4 (HS) 

not measured 

free st. in-plane 

this work LS 

LPCVD 

SiiNi.i 

0.2 GPa 

2.5 

0.8 

free st. in-plane 

this work HS 

LPCVD 

Si3N4 

0.85 GPa 

3 

0.8 

free st. in-plane 


TABLE I. Measured values of thermal conductivity {k in W.m“^.K“^) and specific heat (c in 
J/g.K) of silicon nitride having different stoichiometry and/or different stress. Our results are in 
accordance with most of the studies. LS is for Low Stress, HS for High Stress and free st. for free 
standing membranes. 


3 W.m“^.K“^. These values are in accordance with most of the in plane values of thermal 
conductivity measured which are displayed in table I. Indeed Jain and Goodson^^ have 
measured the in-plane thermal conductivity of 1.5 /rm thick SiN specimens to be about 
5 W.m“^.K“^. At the nanoscale, Sultan et al.^^ reported thermal conductivity of 500 nm thin 
hlms as 3-4 W.m“^.K“^ for a temperature range of 77-325 K. For 180-220 nm thick LS nitride, 
Zink and Heilman^® also observed temperature variation of thermal conductivity ranging 
from 0.07 to 4 W.m“^.K“^ from 3 to 300 K. The cross-plane thermal conductivity measured 
by Lee and CahilF^ for less than 100 nm thickness was in the range of 0.4-0.7 W.m“^.K“^ 
showing severely reduced thermal conductivity, which was ascribed to the interfacial thermal 
resistance. Zhang and Grigoropuolos^® also observed anomalous thickness dependence and 
suggested that micro structural defects may strongly influence thermal conductivity. It is 
important to note that none of the above studies measure thermal conductivity as a function 
of the internal stress. 

In order to verify the coherence of our experimental results, we have extracted the specihc 
heat from the variation of the dcu signal versus the frequency. Generally the specihc heat is 
not expected to vary strongly as a function of stress at room temperature®, and consequently 
it is a good test for the experiment. The results for the four different membranes are shown in 
Fig. 5. The temperature variation of the specihc heat is very similar for the four samples. For 
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both 50 and 100 nm thick membranes we observe that the specihc heat tends to be slightly 
higher for the case of low stress sample. But here again, the differences are insignihcant 
and the specihc heat is very similar for all the thicknesses and stress (low and high). The 
Debye temperatures deduced from the heat capacity measurements vary from 620 to 650 K 
depending on the sample which is a little lower than the commonly accepted value^®. Our 
measurements of thermal conductivity and specihc heat demonstrate that no signihcant 
diherences occur for the thermal transport in high and low stress SiN material because even 
with a stress close to 1 GPa, no modihcation of the phonon thermal conductivity can be 
observed. 



FIG. 5. (colour online) Specific heat measurement of 50 nm and 100 nm thick of both LS and HS 
silicon nitride from low temperature (10 K) to room temperature. 


IV. DISCUSSION 

High stress silicon nitride mechanical devices exhibit remarkable Q factors: inverse quality 
factors Q~^ are two to three orders of magnitude lower than those of amorphous Si 02 from 
4 K up to room temperature^. The true origin of mechanical dissipation in stressed SiN is 
still unknown but could have connections with the thermal properties^’^^’^^. Even though 
amorphous solids are by nature diverse in composition, these materials are characterized 
by a universal behaviour of the thermal conductivity and mechanical dissipation at low 


10 



temperature (between 0.1 and 10 This universal behavior was initially reported 

by Zeller and Pohl^° and described in terms of a phenomenological model which takes into 
account the contribution from defects referred to as two-level systems (TLS)^^’^^. The model 
does reproduce the data, but the universality appears as a surprising coincidence which 
continues to puzzle physicists^’^^. 



FIG. 6. (colour online) Mean free path A of measured samples calculated using experimental data 
of specific heat and thermal conductivity. The dashed line shows the estimation of the mean free 
path using the Debye specific heat. 

In the theoretical work by Wu and Yu®, the starting point is to consider that the stress 
(bond constraints, impurities, local defaults or even external strain) can modify either the 
TLS barrier height V or the coupling between TLS and phonons denoted by 7 . In this model, 
it is predicted that the modification of V and 7 (by taking into account the amplitude of 
the stress in stoichiometric SiN) will have a significant effect on the thermal conductivity 
and mechanical dissipation. Let us discuss the two cases separately. First, when the barrier 
height is affected a difference between the thermal conductivity in low and high stress should 
be seen as the temperature is reduced with a factor close to five at 50 K. This is clearly not 
observed in our measurements since the thermal conductivity of the HS and LS membranes 
are very similar. We can only point out that around 50 K the thermal conductivity is 
slightly different between 50 nm and 100 nm samples, a behaviour that can be attributed to 
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a reduction of mean free path in the thinner membrane. Secondly, in the case of the coupling 
between phonon and TLS (given by the parameter 7 ), an effect even larger is expected with 
a thermal conductivity a factor of ten higher in the HS SiN as compared to the LS at room 
temperature. This could be indeed very interesting for practical applications. Even though 
the stoichiometry is not strictly identical between the low and high stress membranes, we do 
not observe such a big difference in thermal conductivity. This has to be drawn closer to the 
mechanical measurement performed at 4 K, which also did not present any large differences 
in mechanical damping between HS and LS devices. 

We thus demonstrated negligible effect of stress on the thermal conductivity and mechan¬ 
ical dissipation in amorphous SiN. We conclude that the hypothesis of TLS which barrier 
height V or coupling constant 7 is affected by stress does not apply to these materials in 
the present stress range. We also underline that the values of thermal conductivity we have 
measured for both high stress Si 3 N 4 and low stress SiN membranes are in perfect accordance 
with most of the values already published (see table I). In order to highlight the low temper¬ 
ature particularities of the phonon conductivity in these thin membranes, it is particularly 
important to discuss the temperature variation of the mean free path^^. Fig. 6 shows the 
phonon mean free path in the membranes determined from the kinetic equation A = Sk/Cvg, 
Vg being the Debye speed of sound. It has been shown in the past that this equation can 
be used even at room temperature for amorphous materials by Pohl and co-workers^®. At 
300 K all curves (with the exception of 50 nm LS) approach the same limit which is two 
times higher than the inter-atomic spacing (0.25 nm for amorphous SiN). This is in very 
good agreement with previous thermal analysis^’^. As the temperature decreases, the mean 
free path increases rapidly to reach the order of ten nanometers at 20 K. As it can be seen 
in Fig. 6 , it is reasonable to ascribe the difference of thermal transport below 200 K to a 
reduced mean free path in the thinner membranes. 


V. CONCLUSIONS 

The thermal conductivity has been measured on silicon nitride membranes having low 
and high stress. The objective was to search for any effect of internal stress on the phonon 
thermal conductivity and mechanical dissipation. Even though very high stress (of the order 
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of 1 GPa) has been evidenced in suspended stoichiometric SiN membranes by nanomechan¬ 
ical measurements, it has been shown using very sensitive Sea technique that the thermal 
conductivity was not affected. Besides, mechanical dissipation is almost independent of 
stress, even though high Qs are obtained in HS structures in accordance with the ’’dissipa¬ 
tion dilution” model. This rules out a scenario of strong increase of thermal conductivity 
(and concomittantly a strond decrease of mechanical dissipation) with the presence of stress 
proposed recently by Wu and Yu®, either through the increase of the barrier height of two 
level systems or through the decrease of the coupling between TLS and phonons. We also 
show that the thermal properties of the most commonly used silicon nitride materials are 
equivalent. We then express doubts about the possible use of stress in thermal engineering 
in amorphous materials. 
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